Metal-insulator transition caused by the coupling to 
localized charge-frustrated systems under ice-rule local constraint 
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We report the results of our theoretical and numerical study on electronic and transport properties 
of fermion systems with charge frustration. We consider an extended Falicov-Kimball model in 
which itinerant spinless fermions interact repulsively by U with localized particles whose distribution 
satisfies a local constraint under geometrical frustration, the so-called ice rule. Electronic states of 
the itinerant fermions are studied by approximating the statistical average by the arithmetic mean 
over different configurations of localized particles under the constraint. We numerically calculate 
the density of states, optical conductivity, and inverse participation ratio for the models on the 
pyrochlore, checkerboard, and kagome lattices, and discuss the nature of metal-insulator transitions 
at commensurate fillings. The results are compared with exact solutions for the models on Husimi 
cacti as well as with numerical results for completely random distribution of localized particles. 
As a result, we show that the ice-rule local constraint leads to several universal features in the 
electronic structure common to different lattice structures; a charge gap opens at a considerably 
small U compared to the bandwidth, and the energy spectrum approaches a characteristic form in 
the large U limit, that is, the noninteracting tight-binding form in one dimension or the 5-functional 
peak. In the large U region, the itinerant fermions are confined in the macroscopically-degenerate 
ice-rule configurations, which consist of a bunch of one- dimensional loops: We call this insulating 
state the charge ice. On the other hand, transport properties are much affected by the geometry 
and dimensionality of lattices; e.g., the pyrochlore lattice model exhibits a transition from a metallic 
to the charge-ice insulating state by increasing U , while the checkerboard lattice model appears to 
show Anderson localization before opening a gap. Meanwhile, in the kagome lattice case, we do not 
obtain clear evidence of Anderson localization. Our results elucidate the universality and diversity 
of phase transitions to the charge-ice insulator in fully frustrated lattices. 

PACS numbers: 71.10.Fd, 71.27.-fa, 71.30.-fh 



I. INTRODUCTION 

Charge frustration has recently attracted considerable 
attention, since it underlies many intriguing phenom- 
ena such as electronic ferroelectricity, multifcrroicity, and 
novel metallic state in quantum melting of charge or- 
der [H-Q. In these systems, frustration often prevents the 
system from stabilizing a long-range order, and results in 
a ground state with macroscopic degeneracy. The degen- 
erate manifold is highly sensitive to perturbations, such 
as couplings to other degrees of freedom, quantum and 
thermal fluctuations, and external stimuli such as appli- 
cation of pressure or external fields. This is a source of 
fascinating phenomena in charge frustrated systems. 

In the case of localized electrons, the problem is well 
described by the Ising models on frustrated lattices, in 
which spin up and down represent charge rich and poor 
states. A classic example is found in a charge ordering 
on the frustrated pyrochlore lattice, which was argued as 
the origin of the metal-insulator transition in magnetite 
Fe304 [j]. It was pointed out that the problem can be 
mapped onto an antifcrromagnetic Ising model on the py- 
rochlore lattice, and the strong frustration precludes the 
nearest-neighbor Coulomb repulsion to stabilize a charge 
ordering . The degenerate ground-state configurations 
consist of arbitrary network of locally-correlated tetrahe- 
dra, in which every tetrahedron has two charge-rich and 



two charge-poor sites. This local constraint is equivalent 
to the one discussed for proton configurations in water 
ice 0], and is called the ice rule. More recently, a 
magnetic analogue, the spin ice, was discovered in sev- 
eral pyrochlore oxides [1,0, which has been promoting 
the understanding of the ice-rule physics [l0| . 

One of the most striking features of such ice-rule sys- 
tems is cooperative nature of the spatial correlation. It 
was argued that the local constraint brings about a hid- 
den gauge structure: The ice rule corresponds to a zero- 
divergence condition in terms of a notional electric (or 
magnetic) field, which leads to a dipolar correlation in 
the ice- rule variables [ll|, 113 • Systems under the local 
constraint arc not simply disordered but have rather co- 
operative nature. 

Recently, the effect of ice rule on electron itinerancy 
has also attracted increasing interest. It is highly non- 
trivial how the cooperative nature of the ice-rule degen- 
erate manifold affects electronic and transport proper- 
tics. Experimentally, anomalous transport phenomena 
were observed in several pyrochlore-based compounds 
and their relation to the ice-rule degeneracy has been 
discussed. For example, peculiar magneto-transport phe- 
nomena were found in hybrid compounds of itinerant d 
electrons and localized / moments in Nd2Mo207 [l^ and 
Pr2lr207 [3, [l^, in which the / moments potentially 
have spin-ice type correlations. Spin-ice like frustration 
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in charge and orbital degrees of freedom might also be 
relevant to a heavy-mass behavior and a metal-insulator 
transition under pressure in a mixed-valence compound 
LiV204 [T6l420| . Theoretically, although several interest- 
ing aspects were pointed out, such as a fractional charge 
excitation [2lj and a lifting of ice-rule degeneracy by ki- 
netic motion of electrons [l^, [2^, much less is known 
about the role of ice rule in itinerant systems, compared 
to that in localized systems. 

The aim of this study is to elucidate how the cou- 
pling to a localized system with the ice-rule local con- 
straint affects electronic and transport properties of itin- 
erant electrons. We try to reveal the fundamental and 
universal aspects of such systems, instead of explain- 
ing each specific experimental data in real compounds. 
We examine this problem by considering an extended 
Falicov-Kimball model on geometrically frustrated lat- 
tices, such as the pyrochlore, checkerboard, and kagome 
lattices. The model describes itinerant spinless fermions 
interacting with localized particles whose spatial distri- 
bution obeys the ice rule. We numerically calculate elec- 
tronic properties of the model by taking the arithmetic 
mean over different ice-rule configurations of localized 
particles, instead of the statistical average. This is an ap- 
proximation that becomes exact when different ice-rule 
configurations have the identical Boltzmann weight: It is 
indeed the case of the models on Husimi cacti for which 
we obtain exact solutions for comparison. We discuss to 
what extent the approximation holds in the other lattice 
systems. 

Our main result is the clarification of universality 
and diversity of the nature of metal-insulator transition 
caused by the coupling to the ice-rule manifold. Wc find 
that the electronic structure exhibits several universal 
features irrespective of the lattice structures. In partic- 
ular, the systems commonly show a gap opening as in- 
creasing the repulsive interaction at a commensurate fill- 
ing. The gapped insulator in the strong coupling region 
exhibits a peculiar electronic state; itinerant fermions arc 
localized in the specific ice- rule configuration — we call it 
the charge-ice insulator. In contrast, transport properties 
are sensitively dependent on the lattice structure. For the 
three-dimensional pyrochlore-lattice model, our numeri- 
cal results indicate a direct phase transition from a metal 
to the charge-ice insulating state. On the other hand, the 
two-dimensional checkerboard lattice appears to exhibit 
Anderson localization, namely, a metal-to- insulator tran- 
sition without gap opening, before going into the charge- 
ice insulator. For the kagome lattice, however, the results 
do not show any clear sign of Anderson localization. We 
discuss in detail the universal and diverse aspects of the 
effect of ice rule on electron itinerancy through the com- 
parison with the exact solutions for the cactus models 
as well as the numerical results for completely random 
distribution of localized particles. 

The organization of this paper is as follows. In Sec. |lTl 
we introduce models and methods. After introducing the 
model Hamiltonian and the ice-rule limit in Sec. Ill Al 



we describe how we compute the electronic properties in 
Sec. Ill Bl and III Ul The validity of the approximation 
introduced in the calculation is examined in Appendix 
A. We present the various lattice structures in Sec. IIIDI 
and describe the exactly-solvable cases, Husimi cacti, in 
Sec. Ill El In Sec. IIIII we show our results on electronic 
and transport properties of the models on the pyrochlore, 
checkerboard, and kagome lattices, in comparison with 
the exact solutions for Husimi cacti and the numerical 
results for uncorrelated random distributions. The re- 
sults are discussed in detail in Sec.|IVl Sec.|V]is devoted 
to summary. 



II. MODEL AND METHOD 
A. Model 

We start with an extended Falicov-Kimball model [l^] 
(ex-FK model) on the pyrochlore lattice [Fig. [TJa)] , 

(id) 

+ uY,nKn{ ~1) + Vj2n{nl (1) 

where Ci (cj) is an annihilation (creation) operator for 
an itinerant spinless fermion at site i and is the num- 
ber operator, = c|ci. n{ denotes the number of lo- 
calized classical particles at site i; n{ = or 1. The 
sum runs over the nearest-neighbor sites. The first 
term describes the hopping of itinerant fermions, the sec- 
ond term is the repulsive interaction between itinerant 
fermions and localized particles {U > 0), and the third 
term represents the nearest-neighbor interaction between 
localized particles. We take an energy unit as t = 1 ex- 
cept when explicitly shown. 

Hereafter we focus on the "ice-rule limit" : (i) the total 
number of localized particles is fixed at N/2, with N 
being the total number of sites, and (ii) the interaction 
between localized particles is taken to be positive infinity, 
i.e., V/t — ^ oo. In this limit, the localized particles are 
distributed over the system with satisfying the ice-rule 
constraint; two out of four sites are occupied by localized 
particles in every tetrahedron, as exemplified in Fig.[lja). 
We call the model in the ice-rule limit the ice-rule model 
hereafter. 

Later, we extend the model ([T]) to other lattices, the 
checkerboard and kagome lattices in Sec. IIIDI We also 
consider the models on the Husimi cactiin Sec. Ill El for 
which the analytical solutions are available. We will in- 
troduce a similar ice-rule appropriately in each case. 
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B. Arithmetic mean approximation within the 
ice-rule manifold 



For the model ([T]) , the expectation value of an observ- 
able A is given by 



Tr/TrcAexp(~/3i?) 
Tr/Trcexp(-/37?) ' 



(2) 



where Trc (Tr/) is the trace over configurations of itin- 
erant fermions (localized particles) and (3 is the inverse 
temperature. For a given configuration {n{}, the Hamil- 
tonian ([T]) is reduced to a one-body Hamiltonian with 
onsite potential: 



H{{n{}) = -t 



+ h.c.) + U,'i 



(3) 



where the binary potential is given by 



= U{ni --)= + 



U 
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2 



U 
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(4) 



corresponding to 7i{ = 1 or 0. Since we consider the ice- 
rule limit here, the spatial distribution of Ui obeys the 
ice-rule constraint; two sites of every tetrahedron being 
+U/2 and the other two being —U/2. Then, the expec- 
tation value (A) is rewritten as 



(A) 



E 



{n }^ice 



{^g\A\ijg)eXp{-PEg) 



{n . }Gice 



(5) 



where En 



Eg{{nf}) and |V<,) = |V'<,({<})) are the 
energy and many-body eigenfunction of the ground state 
of iJ({n{}); =i?,({nf})|V',({nf})). 
Here, we assume the temperature to be sufficiently low, 
and ignore the contribution from excited states. The 
sum over {n{} £ ice is taken for all configurations of the 
binary onsite potential Ui that obey the ice rule. 

Then, by computing the eigenvalues and eigenstates for 
all the ice-rule configurations, one obtains the statistical 
average of observables. However, it is virtually impossi- 
ble to take the sum in Eq. ([5]) since the number of ice-rule 
configurations grows exponentially with increasing the 
system size 1.5^/ 2) [1 HI). To avoid this difficulty 
and to extract the essential physics of itinerant electrons 
coupled with the ice-rule localized variables, we replace 
the statistical average by the arithmetic mean with omit- 
ting the Boltzmann weight in Eq. ([S])- This corresponds 
to the assumption that Eg{{n{}) does not depend on dif- 
ferent ice-rule configurations of onsite potential Ui, i.e., 
Eg{{nj}) ~ Eq. Moreover, we calculate the average by 
sampling the ice-rule configurations randomly. Namely, 



we approximate Eq. ([5|) by 



(A) 



E{„/}6iceexp(-/3i^o) 



-'''ice ■ T 1 
2^1 m— 1 



— E E<"^i^i"^)' 



(6) 
(7) 
(8) 



samp 



= 1 m=l 



where A^ico ^ 1.5^/^ is the total number of ice- rule con- 
figurations of localized particles, A'samp is a number of 
samples in actual calculations, and Nc is the total num- 
ber of itinerant fermions. Here, |m) are the one-particle 
eigenstates of H{{n{}), and are sorted in ascending order 
of the eigcnenergies. 

The arithmetic mean in Eq. ([7]) becomes exact if the 
Boltzmann weight in Eq. ([5]) is identical for all differ- 
ent ice-rule configurations of localized particles in the 
ground state. This is indeed the case for models on 
Husimi cacti [l^l as discussed in Sec. IH El In these mod- 
els, because of the loopless structure of lattices, the ice- 
rule configurations are all topologically equivalent, which 
results in the identical Boltzmann weight. For the py- 
rochlore lattice case, however, the Boltzmann weights 
arc not identical for different configurations; neverthe- 
less, the differences turn out to be very small because of 
the structure of degenerate ice-rule manifold. In fact, we 
find that the difference of the total energy is typically in 
the order of 10~H in the entire region of U/t. Detailed 
discussions will be given in Appendix [X] This suggests 
that the approximation in Eq. ([7]) gives quantitatively 
reasonable results in the situation that we are interested 
in, where the ice-rule manifold is well preserved and the 
system takes all the ice-rule states without selecting a 
unique ground state or a submanifold if any. In other 
words, the arithmetic mean provides a tractable and rea- 
sonable tool to investigate the effect of the ice-rule de- 
generacy on the electronic state of itinerant fermions. 

Through the procedure above, the many-body prob- 
lem in Eq. ([T]) is reduced to a one-body problem with 
onsite ice-rule potential in Eq. ^ . For a reference to the 
ice-rule model, we also consider a one-body model with 
random potential. In this case, we consider completely 
random configurations of N/2 localized particles without 
any spatial correlation in the Hamiltonian ([3]) . The com- 
parative study illuminates the effect of local correlation 
in the ice-rule model, as we will see later. 

In general, a special configuration of the potential may 
lead to a characteristic electronic state. For example, it 
was pointed out that, even in the one-dimensional 
locally-correlated potential can drive the system delocal- 
ized |27[. In the following, we elucidate the effect of local 
correlation brought by the ice rule in higher dimensions, 
in comparison with the random cases without any spatial 
correlation. 
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C. Numerical calculations 

We calculate the electronic properties of model ([T]) nu- 
merically by using Eq. In the calculation, we need 
to generate iVsamp samples of ice-rule potential configu- 
rations and obtain the eigenstates for each sample. To 
generate different ice-rule samples sequentially, we em- 
ploy the so-called loop algorithm [H, [2^ . To retain the 
statistical independence, we applied 10^ loop updates be- 
tween the samples. This is much larger than the "auto- 
correlation time" T, typically in the order of tens [soj . 
On the other hand, we also examined the convergence of 
Eq. ([5]) as to Nsamp- We checked iVgamp dependence and 
conclude that average over A^^^amp = 40 gives converged 
results with enough precision for the following discus- 
sions. Statistical errors are fairly small for the density 
of states (DOS), which are comparable to the width of 
curves in the plots below. For other quantities, the errors 
are explicitly shown in the plots. Once the sample set is 
obtained, we calculate the eigenenergies and eigenstates 
by the exact diagonalization of the Hamiltonian given by 
Eq. ©. 

In the following sections, we discuss the electronic and 
transport properties by calculating DOS and the opti- 
cal conductivity. DOS is obtained from the eigenener- 
gies by taking the histogram over A^samp samples with 
an energy width Ae = 0.02. The energy gap is directly 
calculated by the arithmetic mean of the energy differ- 
ences between the highest occupied level and the lowest 
unoccupied level. The optical conductivity is calculated 
by the standard Kubo formula: 



(9) 

where /(e) is the Fermi distribution function and 

= ("a" ■ ^],k){clcj - c]cfe) (10) 

is a current operator in the /i direction is assigned 
for each case below), which is constructed in a standard 
way from a polarization operator in order to satisfy the 
continuity equation j3l| . Here, is the unit operator in 
the ^ direction, t is the transfer integral for the nearest- 
neighbor sites, and dj^k is the geometrical vector from jth 
to fcth site. The sum is taken for all the nearest-neighbor 
pairs. To calculate the optical conductivity at T — > 0, 
we set /3 = 10^. As a measure of the metallicity, we also 
calculate the low-energy weight of cr(a;) defined by 

w ^ <t{uj) duj (11) 
Jo 

with ujQ being a cut off. 



To further examine the metallicity of the system, in 
particular, to detect Anderson localization, we also com- 
pute the inverse participation ratio (IPR) [32|. IPR is 
defined by 

p-l^5]l^(R0l^ (12) 

I 

where </>(R;) denotes an eigenfunction at site R;. IPR ex- 
trapolated to the infinite system size N gives a measure of 
the localization of the eigenfunction: If (f) is an extended 
state, IPR behaves as oc when N 00, 

while if (f) is localized, converges to a nonzero value 
as iV — >■ 00. Usually, IPR is calculated for the state right 
at the Fermi level, i.e., for the highest occupied state. In 
the following calculations, however, we focus on a special 
filling at which DOS shows divergence at the Fermi level 
due to the degeneracy of a huge number of states. (See, 
e.g.. Fig. [31) In such cases, IPR right at the Fermi level 
is not well defined because of the degeneracy. Hence we 
measure the metallicity by calculating IPR for the states 
just above the flat bands. The evaluation is conducted by 
the harmonic average of IPR for the states within a small 
energy window with the width of Ae = 0.02, so that the 
averaged IPR goes to zero as TV — )■ 00 if the states in the 
energy window include at least one extended state [s^ . 

D. Lattices 

The model ^ is introduced for the pyrochlore lattice 
shown in Fig. [IJa). In the following, we also consider 
similar models on the checkerboard and kagome lattices, 
as shown in Figs.[IJb) and[TI^c). These lattice structures 
are two-dimensional cousins of the three-dimensional py- 
rochlore lattice; the checkerboard lattice is a (001) pro- 
jection of the pyrochlore lattice, while the kagome is a 
(111) plane. The checkerboard lattice shares its geo- 
metrical imit, a tetrahedron, with the pyrochlore lattice, 
while the unit of the kagome lattice is a triangle. On the 
other hand, the kagome lattice shares a global geometri- 
cal feature with the pyrochlore lattice, i.e., the smallest 
loop composed of the geometrical units is a hexagon as in 
the pyrochlore case; whereas it is a square in the checker- 
board case. Through the comparative study among these 
different lattices, we examine the effect of differences in 
dimensionality, geometrical unit, and their connection. 

For the checkerboard lattice, the ice-rule local con- 
straint is applied in the same manner as for the py- 
rochlore lattice. Namely, we consider N/2 localized par- 
ticles, and place two of them on each tetrahedron, as in 
Fig. [Ijb). On the other hand, for the kagome lattice, we 
consider an ice-rule type constraint by distributing N/3 
localized particles so that one site is occupied and the 
other two are unoccupied in each triangle, as exempli- 
fied in Fig. [TJc). We call this rule the "kagome ice rule" 
hereafter. 

All three lattice structures share common features; 
they consist of corner sharing network of geometrically- 
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frustrated units, tetrahedra or triangles. These units are 
called the complete graph, i.e., a graph in which all ver- 
tices arc connected with all the other vertices. To illu- 
minate the role of the common features, we also consider 
a set of variants for these lattices, that is, the so-called 
Husimi cacti of the complete graphs, as discussed sepa- 
rately in the next section. 



E. Exactly solvable models on the Husimi cacti 

In addition to the lattice structures introduced in the 
previous sections, we also consider modified structures 
composed of the same geometrical units — the Husimi 
cacti of tetrahedra and triangles. They arc the analogues 
of the Bethe lattice composed of corner-sharing tetrahe- 
dra or triangles, which we call here the tetrahedron or 
triangle Husimi cactus, respectively [Figs.[Ud) and[Ue)]. 
The Husimi cacti share two important structural features 
with the original lattices; the geometrical units and their 
corner-sharing network. A difference is in the global con- 
nection of the units: As mentioned in Sec. IHDl all the 
three lattices have loops running across different geomet- 
rical units, but the Husimi cacti do not have such global 
loops. Thus, the Husimi cacti are the loopless variants 
of the original lattices. 

A crucial advantage of considering the Husimi cacti 
is that the ex-FK model in the ice-rulc limit is exactly 
solvable. The benefit is a direct consequence of the 
loopless structure, which makes different ice-rule con- 
figurations topologically equivalent. The derivation of 
the exact solution for the tetrahedron Husimi cactus 
was already reported in Ref. [l^l, and it is straightfor- 
ward to extend the method to the triangle Husimi cactus 
case. Hence, here we only show the final results. For 
the cacti models, the local Green's function at site i, 
Gi{e) = t{i\[e — 'H{{n{}) + iS]~^\i), depends not explic- 
itly on i but only on the value of Uf, i.e., it is possi- 
ble to write Gi = G± corresponding to the sites with 
Ui = ±C//2. (We include t in the definition of G to sim- 
plify the following expressions.) G± are obtained as 



1 



G± = — -T(eT 

9± t 2 ' 



(13) 



where g± are the solutions of the following recursive equa- 
tions for each case: (i) for the tetrahedron Husimi cactus 
with Y.^ n{ - N/2, 

g±(l-.9T) + 2gT(l-g±) _ 1 .14^ 

l+5^(l-25±) t^'^ 2> g±' ^ ^ 

and (ii) for the triangle Husimi cactus with ^ - n{ = N/3, 



2g- 



1 



1 



5+.9--1 2^ 5 



1 



(15) 
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FIG. 1. (color online). Lattice structures of (a) pyrochlore, 
(b) checkerboard, (c) kagome lattices, and (d) tetrahedron 
and (e) triangular Husimi cacti. DOS for the noninteracting 
tight-binding model on each lattice structure are shown in 
(f)-(j). In the lattice figures, filled (open) circles denote the 
sites unoccupied (occupied) by localized particles. The con- 
figurations represent examples of the ice-rule configurations. 
Loops connecting the unoccupied sites are shown by the bold 
lines. See text for the details. 



Equation with the solutions of Eqs. ([H]) and 
gives the exact local Green's functions of the ex-FK mod- 
els given by Eq. ([T]) on the Husimi cacti in the ice-rule 
limit [U, [iH. As we will show in the following sections, 
the solutions for Husimi cacti models give good references 
to the original lattice models. 
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III. RESULTS 

In this section, we discuss the effect of ice-rule con- 
straint on electronic and transport properties of models 
on different lattice structures one by one. Comparisons 
with the results for models with random potential and 
the exact solutions for the Husimi cacti are also given. 

A. Pyrochlore lattice 

Figure [5] shows DOS for the pyrochlore lattice case. 
The left column (A1)-(A5) shows the results for the 
ice-rule models, the middle (B1)-(B5) for the random- 
potential, and the right (C1)-(C5) for the exact results 
for the tetrahedron Husimi cactus. Different rows cor- 
respond to the data at different values of U. The bold 
(dotted) curves represent the site-resolved DOS, p+ 
for Ui = +U/2 {~U/2) sites, and the thin curves denote 
the total DOS p ^ p+ + p-. The numerical results for 
the pyrochlore lattice (left and middle columns) are ob- 
tained for 3^ superlattices of 4 x 8^ sites. The results for 
the tetrahedron Husimi cactus (right column) are calcu- 
lated from Eqs. (HH) and ([Ti|) . 

In the pyrochlore case, at [/ = 0, the energy levels con- 
sist of two flat bands and two dispersive bands [Fig.[T](f )]. 
The flat bands give the (5- functional peak at e = 2. Mean- 
while, the dispersive bands are equivalent to those of the 
diamond lattice. This equivalence can be understood 
as a result of line-graph correspondence. The dispersive 
bands form a continuum spectrum for —6 < s < 2 with 
two semimetallic dips; one is at e = — 2, where two dis- 
persive bands touch with each other, and the other is at 
e — 2, where the higher dispersive band touches the flat 
bands. At half-filling ~ -^72; the Fermi energy is 

located at e = 2, right at the latter semimetallic point; 
i.e., two dispersive bands are fully occupied and the fiat 
bands are empty. 

By switching on U, the flat bands are perturbed to be 
broadened, resulting in a spectrum sandwiched by two 
divergences in DOS, as shown in Figs.^Al) andjU^Bl). 
The Fermi level at half filling is pinned at the lower- 
edge divergence for both the ice-rule and random cases. 
It is worthy to note that, in the ice-rule case, a cusp- 
like structure appears between the two divergences, as 
indicated by arrows in Figs. [51JA1)-(A4). We return to 
this point below. 

As U increases further, an energy gap starts to open 
both for the ice-rule and random cases. However, the 
critical values of gap opening, Uc, are largely different 
between these two cases. In the results for icc-rule case, 
there is a clear gap with U >3 [Figs. [U A3) and[i;A4)], 
whereas there still remains a small DOS at the Fermi level 
at half filling in the random case until U ^ 5 [Figs.[21[B3) 
and [2l^B4)] . The behavior of the gap opening near Uc is 
shown for the ice-rule model in Fig. [3] U dependence of 
the energy gap at half filling is summarized in Fig. ^a). 
For the random case, a gap appears at [/ ~ 6, i.e., when 



the potential value becomes comparable to the band- 
width. On the other hand, the critical value of U in 
the ice-rule case is estimated as Uc ~ 2.3(2). Such a 
small Uc compared to the bandwidth is considered to be 
characteristic of the correlated ice-rule configurations. 

In the ice- rule case for U > Uc, the system is insulating 
at half filling, where itinerant fermions are excluded from 
the sites with Ui = +U/2 and localized in the ice-rule 
configurations at Ui = —U/2 sites. We call this insu- 
lating state the charge-ice insulator (2^. In the large U 
limit, DOS for the charge-ice insulator approaches a pair 
of DOS for the one-dimensional (ID) tight-binding model 
centered at e = ±U/2, as shown in Fig.[2|A5). This is in 
sharp contrast to the asymmetric featureless DOS for the 
random case in Fig. [2IB5). The upper-edge divergences 
of the two ID-like bands at e ~ ±[//2 -I- 2 come from 
the two divergences in the perturbed flat bands in the 
small U region described above. On the other hand, the 
lowcr-cdge divergence in the upper band at e ~ U/2 — 2 
develops from the cusp-like feature seen already at e ~ 2 
in the small U region. The characteristic ID-like fea- 
ture is explained by the fact that, in the ice-rule case, 
the sites with Ui ~ +U/2 and —U/2 form ID loops of 
equivalent onsite potential [see Fig. [TJa)]. The ID loops 
have many different lengths from the shortest six site 
to infinite length; the average over the lengths for dif- 
ferent ice-rule configurations results in the ID- like DOS. 
The form of DOS shows slight deviations from that for 
the ID tight-binding model; a dicernible feature is the 
spikes in the spectrum, which is presumably due to the 
finite length of loops. Thus, the gapped insulating state 
at large U is the charge-ice insulator in which itinerant 
fermions are confined in the ID loops. 

The interesting evolution of DOS with increasing U for 
the ice-rule case is well reproduced by the tetrahedron 
Husimi cactus model [Figs. [2i;Cl)-(C5)] H. This model 
shows a transition to the charge-ice insulator at Uc ~ 
2, much smaller than the bandwidth, similarly to the 
pyrochlore lattice model. Furthermore, in the large U 
limit, DOS converges to exactly the same form for ID 
tight-binding model centered at ±C//2. This comes from 
the fact that the cactus lacks finite-length loops which are 
present in the pyrochlore lattice [Fig. [ijd)]. The cactus 
also well reproduces other important features, such as 
the overall form of DOS at U ~ (dispersive part plus 
flat bands [23|), [Fig. [Tfi)], the split of the flat bands by 
U, and the cusp-like feature. The agreement indicates 
that these universal features are owing to the peculiar 
corner-sharing geometry of tetrahedral units. 

Now we turn to the transport properties. Figure [J] 
shows the optical conductivity at half filling calculated 
along the [111] direction for 4^^ superlattices of 4x4'^ sites. 
The solid (dotted) curves show the results for the ice-rule 
(random) case. As U increases, cr(w) evolves very differ- 
ently for the ice-rule and random cases. In the random 
case, (t(w) changes slowly with developing a dip at cj 0. 
In contrast, in the ice-rule case, the change is more dras- 
tic. First of all, a gap opens for U > Uc, consistent with 
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FIG. 2. (color online). DOS of itinerant fermions for the pyrochlore lattice models with ice-rule constraint (A1)-(A5) and with 
the random potential (B1)-(B5); DOS for the tetrahedron Husimi cactus model with the ice-rule constraint (C1)-(C5). The 
insets of (B2), (B3), and (B4) show the enlarged figures of the main panels in the vicinity of the Fermi level at half filling. Bold 
(dotted) curves represent the partial DOS at the sites with potential +U/2 (— C//2). Thin curves represent the total DOS. The 
Fermi level for half-filling case is indicated by the vertical dashed lines. Error bars are within the width of the curves. Arrows 
in the ice-rule cases indicate the cusp-like structures discussed in the text. 



the DOS in Figs. [H; A3) and[li;A4). Furthermore, the 
low-energy part of (t{uj) does not decrease in a mono- 
tonic way before the gap opens; Comparing Figs. |4l^a) 
and mjb) , we observe an increase of the low-energy part 
as [/ Uc- 

This enhancement is more clearly seen in the low- 
energy weight w defined in Eq. The results are 
plotted in Fig. IHJb). Here we take the cutoff luq = 0.095; 
the results are qualitatively independent of when it 
is small enough. The data clearly indicate the increase 
of Ticff with U for the ice-rule case, showing a maximum 
at J7 ^ 2, and sharply drops at U ~ C/c, in contrast 
to the featureless gradual decrease in the random case. 
This characteristic behavior might be due to the pecu- 
liar semimetallic behavior at J7 = 0; the metallicity is 
suppressed in the small U region since there are less low 



energy states available due to the semimetallic dip. We 
will discuss this behavior in comparison with the result 
for the checkerboard lattice case in the next section. 

Another peculiar difference in a{uj) between the ice- 
rule and random cases is seen in the large U region, that 
is, a sharp spike only existing in the random case. In 
Fig. md), the data at t/ = 5 exhibit a spike at e ~ 5 
for the random case: The spike comes from the transi- 
tion between the divergences of DOS at e = ±[//2 + 2. 
However, the spike is suppressed in the ice-rule case. 

The suppression is a consequence of quantum phase 
interference among the eigenstates at the divergence of 
DOS. The divergences come from a particular set of 
states with ID character, which are originally included in 
the flat bands at [/ = 0. In the ice-rule case, each wave 
function is confined in a ID potential loop, and has a uni- 
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FIG. 3. (color online). DOS for the ice-rule model on the 
pyrochlore lattice in the vicinity of the Fermi level at half 
filling. 
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FIG. 4. (color online). Optical conductivity calculated by 
Eq. ([9]) for the pyrochlore lattice models at half-filling. Solid 
(dotted) curves represent the results for the ice-rule (random) 
case. Typical errorbars are shown in each figure. 
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FIG. 5. (color onhne). IPR obtained by Eq. ^ for the 
half-filled state in the pyrochlore lattice models with (a) the 
ice rule and (b) random potential. The lines show the extrap- 
olation to — >■ oo. See the texts for details. 



suggests that the system remains metallic until U = Uc 
and the localization occurs simultaneously with the gap 
opening. On the other hand, the extrapolated value in 
the random case becomes nonzero even at U =1. This 
result indicates that the system is insulating because of 
the Anderson localization [3^ . since DOS remains finite 
at the Fermi level [Figs.[31[Bl)-(B4)] and the gap does not 
open until [/ ~ 6 [Fig. [BJa)]. The extrapolated values to 
N ^ (X are summarized in Fig. |6l^c) |34{ . 

In summary, the above results reveal peculiar features 
in the ice-rule model at half filling, in sharp contrast to 
the random model: (i) The charge gap opens at a small 
Uc compared to the bandwidth, (ii) DOS shows a ID form 
in the large-C/ charge-ice insulator, (iii) 17(0;) exhibits a 
rapid change for U in the low-energy part and a suppres- 
sion of a sharp spike in the large U region, (iv) w shows 
a characteristic nonmonotonic behavior as a function of 
U, and (v) IPR suggests that a metal-insulator transition 
occurs at U = Uc, where the gap opens. All these fea- 
tures are characteristic of the locally-correlated ice-rule 
potentials, which are not seen in the random systems. 



form amplitude and an alternating sign as H \ \ • • 

along the loop. According to Eq. (|9]), the contribution to 
(7(0;) from the transition between the two divergences is 
proportional to j , where is an 

eigenstate in the lower (higher) divergence. Since the 
current operator J has a matrix element only between 
neighboring sites, the summation counts the matrix ele- 
ments for a pair (li, hj) for which the corresponding ID 
loops have neighboring sites. There are, however, cancel- 
lations between the matrix elements due to the alternat- 
ing sign of the wave functions. This interference is the 
reason why the sharp spike is suppressed in the ice-rule 
case. 

Difference in the nature of metal-insulator transition 
is also examined by using the harmonic average of IPR, 
defined in the previous section. Figure [S] shows the re- 
sults. In the ice-rule case, IPR is extrapolated to zero as 
iV — > 00 for [/ < Uc, and becomes finite for U > Uc- This 



B. Checkerboard lattice 

Figure [7] shows DOS for the checkerboard lattice mod- 
els. The left column (A1)-(A5) shows the numerical re- 
sults for the ice-rule model, while the right column (Bl)- 
(B5) for the random configuration. The symbols are com- 
mon to those in Fig. [21 The results are obtained for 3^ 
superlattiees of 4 x 28^ sites. 

DOS for the checkerboard lattice at ?7 = consists of a 
dispersive band and a flat band, similar to the pyrochlore 
lattice [Fig. [TJg)]. The dispersive band is equivalent to 
that of noninteracting tight-binding model on the square 
lattice with an energy shift of —2, and the flat band gives 
(5- functional peak at e = 2. The Fermi level at half filling 
is located just below the flat band; the dispersive band 
is fully occupied while the flat band remains empty. The 
system is metallic, in contrast to the pyrochlore case in 
which the system is semimetallic. 
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FIG. 6. (color online). U dependence of (a) the energy 
gap, (b) the low-energy weight of <j{lj), w, and (c) IPR for 
the pyrochlore lattice models at half filling, w is calculated 
by Eq. ([TTl) with taking the cutoff ojo = 0.095, and IPR is 
estimated by the extrapolation to — >■ oo in Fig.O The inset 
of (a) shows an enlarged view in the vicinity of Uc = 2.3(2). 



Despite of the difference in the dimensionality of lattice 
structures, the evolution of DOS with increasing U shows 
many common features with the pyrochlore lattice. The 
flat band is perturbed to be broadened by U and splits 
into two divergences, as shown in Figs.[71[Al) and[71jBl). 
A cusp structure appears in between the two divergences 
for the ice-rule case as indicated by arrows in Figs.[71JAl)- 
(A4). The Fermi level at half filling is pinned at the lower 
divergence. The energy gap opens at a smaller U in the 
ice-rule case than in the random case [Figs. El^AS), (A4), 
(B3), and (B4); Fig. [i Fig. HHa)]: The critical value of 
U is estimated as Uc = 2.7(2), which is markedly smaller 
than the bandwidth [see also Fig.fTlTa)]. DOS at f7 — > cx) 
approaches a similar form to ID tight-binding model, as 
shown in Fig. [71[A5). All these features are commonly 
seen in the pyrochlore lattice model as well as in the 
tetrahedron Husimi cactus in Sec. IIII Al These features 
illuminate the effect of the local ice-rule configuration 
characteristic of the corner-sharing tetrahedra, working 
irrespective of global features such as lattice structures 
and dimensions of the system. 

Similar features to the pyrochlore lattice are also ob- 
served in the optical conductivity. Figure |9] shows the 
optical conductivity for the checkerboard lattice at half 
filling measured in the x direction. The symbols are com- 
mon to those in Fig. 21 and the calculations were done on 
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FIG. 7. (color online). DOS of itinerant fermions for the 
checkerboard lattice models with the ice-rule constraint (Al)- 
(A5) and with the random potential (B1)-(B5). The insets of 
(B2), (B3), and (B4) show the enlarged figures of the main 
panels in the vicinity of the Fermi level at half filling. The 
symbols are common to those used in Fig. [2] 



3^ superlattices of 4 x 10'^ sites. The overall evolution of 
a{uj) with increasing U resembles that of the pyrochlore 
models; In the random dip slowly develops at 

a; and it remains as a pseudogap with increasing U, 
while a gap opens in the ice-rule case corresponding to 
the gap opening of DOS. 

However, when one carefully looks at the low-energy 
part of cr(i^), there is a difference compared to the py- 
rochlore case. Figure [TTJb) shows the low-energy weight 
for the checkerboard lattice models. W for the ice-rule 
model shows a monotonic decrease with increasing U, 
which is qualitatively different from the nonmonotonic 



10 



0.15 



O 

Q 



0.10 - }j 



0.05 - 



0.00 




n r 

l|J=2.67 — 
U=3.00 — 
U=3.33 




J -*''\ I L 



0.2 0.4 0.6 0.8 1.0 1.2 
£ 



FIG. 8. (color online). DOS for the ice-rule model on the 
checkerboard lattice in the vicinity of the Fermi level at half 
filling. 
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FIG. 9. (color online). Optical conductivity for the checker- 
board lattice models at half filling. The symbols are the same 
as in Fig. |4] 



behavior seen in the pyrochlore case [Fig. El^b)]. This 
is presumed to be owing to the difference in DOS at 
U = Q. As is described above, in the pyrochlore case, 
a scmimctallic gap exists at the Fermi level of half filling, 
on the other hand, the checkerboard lattice is metallic 
with finite DOS at the Fermi level. Hence, at J7 = 0, 
(j{uj) is zero in the former case, but divergent in the lat- 
ter. The different behavior of (j{ui) in the small U region 
presumably comes from this difference. 

Further difference between the pyrochlore and checker- 
board models is seen in IPR. Figure [TOl shows IPR for the 
ice-rule [Fig. [TUl^a)] and random case [Fig. fTUTb)] with 
varying U . Unlike the pyrochlore case, the extrapolated 
values of IPR for both ice-rule and random cases appear 
to remain finite at /7 = 2, which is definitely smaller than 
Uc = 2.7(2). This suggests that the wave functions near 
the Fermi level are localized in spite that DOS does not 
show a gap; namely, the system is Anderson insulator for 
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FIG. 10. (color online). IPR for the half-filled state in the 
checkerboard lattice models with (a) the ice-rule and (b) ran- 
dom potential. 
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FIG. 11. (color online). U dependence of (a) the energy gap, 
(b) the low-energy weight w, and (c) IPR for the checkerboard 
lattice models at half filling, w is calculated by Eq. (|lip with 
taking the cutoff uia = 0.095, and IPR is the extraporated 
value to TV — !• oo. The inset of (a) shows an enlarged view in 
the vicinity of Uc = 2.7(2). 



both ice-rule and random cases. This is in sharp contrast 
to the pyrochlore case, where the ice-rule model appears 
to remain metallic until the gap opening at U = Uc- 

To summarize, the checkerboard lattice show several 
features common to the pyrochlore lattice, in particular, 
in the evolution of DOS with increasing U . On the other 
hand, transport properties exhibit differences, suggesting 
the qualitatively different nature of the metal-insulator 
transition in the ice-rule models between the pyrochlore 
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and checkerboard cases. These similarities and differ- 
ences will be discussed further in Sec. IIVI 



C. Kagome lattice 

Next, we move to the kagome-ice model introduced in 
Sec. IIIDI In this model, we focus on the 2/3-filling case 
of itinerant fermions instead of half filling, since we are 
interested in the transition to charge-ice insulating state 
which is expected to occur at 2/3 filling by confinement 
of fermions to — C//2 sites. We also make a comparison 
with the random model, in which — U/2 (4-J7/2) potential 
sites are randomly distributed with keeping the ratio of 
the total number of sites at 2 : 1 . 

Figure [12] shows DOS for the kagome models. Each 
column shows the numerical results for the kagome ice- 
rule case, the corresponding random case, and the exact 
solution for the triangle Husimi cactus model, respec- 
tively from left to right. Each row shows the results for 
varying U . The calculations for the kagome models are 
conducted for 3^ supcrlattices of 3 x 36^ sites. The re- 
sults for the triangle Husimi cactus model are obtained 
by Eqs. ^ and p5|. 

DOS for the kagome lattice model at ?7 = consists of 
two dispersive bands which form the continuum spectrum 
at -4 < e < 2 and a flat band at e = 2 [Fig. ^\\)]. The 
dispersive bands are equivalent to those of honeycomb 
lattice, and a semimetallic dip exists at e ~ —1, where 
the lower two bands linearly cross by forming the Dirac 
points. At 2/3 filling, the Fermi level is located just below 
the flat band and the system is metallic similar to the 
checkerboard lattice model. 

With switching on J7, the flat band is broadened and 
the lower edge shows a divergence: The Fermi level at 
2/3 filling is located at the divergence [Figs. [rsT AD and 
[TST BI)]. With further increasing [/, in the kagome ice- 
rule case, the gap starts to open at a fairly small Uc = 
1.4(2) compared to the bandwidth [Fig. [l3i;A2); Fig.[13l 
Fig.fTBTa)]. On the contrary, for the random case, the gap 
does not open until C/ ^ 4 — 5, while a quasigap feature 
develops forU >3 [Figs.[ni;B2)-(B4); Fig.IUJa)]. These 
behaviors are qualitatively similar to the pyrochlore and 
checherboard lattice models. 

On the other hand, the evolution of DOS toward 
large U limit is different from the preceding two mod- 
els due to difference in the geometrical unit of the lat- 
tice. In the kagome ice-rule model, the bandwidth of 
the upper band shrinks with increasing U for U > Uc 
[Figs. [T2JA2)-(A4)], and finally, becomes a (5-functional 
peak, as shown in Fig. [T2FA5). Correspondingly, a cusp- 
like structure appearing between two divergences is ob- 
scure [Figs. [T^Al)-fA4)]. in contrast to the pyrochlore 
and checkerboard cases. On the other hand, the lower 
band approaches a ID-like form. These behaviors are 
explained by considering the strong coupling limit: In 
contrast to the pyrochlore and checkerboard cases, the 
sites with Ui = +U/2 are disconnected from each other 



in the limit of ?7 — > cx) in the kagome-ice model, while the 
—U/2 sites form ID loops [see Fig. [T]^c)]. Consequently, 
localized states at each isolated +U /2 sites contribute to 
the (5-functional divergence, and the states for ID loops 
with -U/2 potential give the ID-like DOS. 

The evolution of DOS for the kagome-ice rule case is 
well captured by the triangle Husimi cactus model, as 
shown in Figs. [T2jCl)-(C5). In this cactus, the critical 
value of U for the gap opening is Uc ^ I, which is a half 
of Uc = 2 for the tetrahedron Husimi cactus. 

Another peculiar feature of DOS for the kagome model 
is the evolution of characteristic structures in the disper- 
sive bands, such as the semimetallic dip between the two 
dispersive bands (at £ = — 1 for C/ = 0), and the van-Hove 
singularity (at e = — 2 for [/ = 0). In the kagome ice-rule 
case, these features remain to be clearly seen even for 
U = 5, while they are smeared out already for U = 3 in 
the random case. Note that similar tendency can be seen 
in the pyrochlore and checkerboard cases. 

Next we examine transport properties of this system. 
The optical conductivity at 2 /3 filling behaves in a qual- 
itatively similar manner to that for the checkerboard 
model, monotonically decreasing with increasing U as 
shown in Fig.[T31 The calculation is done for 4^ supcrlat- 
tices of 3 X 12^ sites. A gap in the low-o; region rapidly 
develops in the kagome-ice case corresponding to the gap 
opening in DOS, whereas ct(oj) changes rather slowly in 
the random case. The low-energy weight w is also shown 
in Fig. I16r b) . w decreases monotonically with increasing 
U and goes to zero at J7 ^ Uc in the kagome ice-rule case, 
similarly to the checkerboard model in Fig. [TTT b'). 

The result of IPR for this system is shown in Fig. [131 
With sufficiently small U, the results for the kagome- 
ice model appear to show the wave functions to remain 
extended, while they appear to be localized with very 
small U in the random case. The extrapolated values of 
IPR shown in Fig. \lE[ c) indicates this difference more 
clearly: The extrapolated values of IPR remain to be 
almost zero for U < Uc in the ice-rule case, while they 
are finite even at C/ = 1 in the random case. These 
behaviors will be discussed further in comparison with 
the pyrochlore and checkerboard cases in Sec. IIVI 

As a consequence, our calculations for the kagome-ice 
model suggest a transition from metal to gapped charge- 
ice insulator without any clear indication of Anderson 
localization. On the other hand, in the random case, the 
Anderson insulator exists before the gapped insulating 
state. The results will be discussed in the next section, 
in comparison with those for the other lattice models. 



IV. DISCUSSIONS 

All the results in Sec. IIIII clearly elucidate peculiar 
electronic properties under the influence of ice-rule cor- 
relation, which are distinctive from those in the com- 
pletely random case. Furthermore, the comparative stud- 
ies among different lattice structures revealed common 
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FIG. 12. (color online). DOS of itinerant fermions for the kagome models with the kagome ice-rule constraint (A1)-(A5) and 
with the random potential (B1)-(B5); DOS for the triangle Husimi cactus model with the kagome ice-rule constraint (C1)-(C5). 
The insets of (B2), (B3), and (B4) show the enlarged figures of the main panels in the vicinity of the Fermi level at 2/3 filling. 
The symbols are common to those in Figs. [2]and[7l 




FIG. 13. (color online). DOS for the ice-rule model on the 
kagome lattice in the vicinity of the Fermi level at 2/3 filling. 



FIG. 14. (color online). Optical conductivity for the kagome 
models at 2/3 filling. The symbols are the same as in Figs. U 
andU 
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FIG. 15. (color online). IPR for the 2/3-filled state in the 
kagome models with (a) the kagome-ice rule and (b) the ran- 
dom cases. 
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FIG. 16. (color online). U dependence of (a) the energy gap, 
(b) the low-energy weight w, and (c) IPR for the kagome 
models at 2/3 filling, w is calculated by Eq. (|11|) with taking 
the cutoff ujo = 0.095, and IPR is the extraporated value to 
— > oo. The inset of (a) shows an enlarged view in the 
vicinity of Uc = 1-4(2) 



features and dilFcrcnccs. We discuss the universality and 
diversity in the follovifing. 

First, we discuss several universal features of DOS 
common to the three models we have investigated 
through this paper. All the models exhibit a consider- 
ably small critical value of U for gap opening, Uc- The 
values of Uc are estimated as 2.3(2), 2.7(2), and 1.4(2) 
for the pyrochlore, checkerboard, and kagome models, re- 
spectively. These values are markedly smaller than the 
bandwidth at U = 0, W; = 8 for the pyrochlore and 



checkerboard lattices, and W = 6 for the kagome lattice. 
On the other hand, for the random case, the energy gap 
evolves around U ^ W for all the three models, reflecting 
the competition between the local binary potential and 
kinetic energy. This contrastive behavior of Uc indicates 
that the ice-rule short-range correlation helps to develop 
a gapped insulating state, compared to the uncorrelated 
random case. 

Second, DOS shows a common feature in the large U 
regime. For all lattice structures, at J7 — > oo, itiner- 
ant fermions are excluded from Ui = +U/2 sites and 
localized at Ui = —U/2 site forming ice-rule charge dis- 
tribution complementary to that of localized particles. 
We call this insulating state the "charge ice" insulator. 
The distribution of fermions is a set of ID loops made of 
-~U/2 potential sites. Reflecting the specific configura- 
tions, DOS approaches a universal symmetric form given 
by the ID tight-binding form or the S function. This be- 
havior of DOS in the large U limit is also qualitatively 
different from those for random models: In the random 
cases, DOS becomes highly asymmetric and broadened, 
irrespective of the lattice structure. 

In addition, we note that the characteristic features in 
the noninteracting band structure are well preserved for 
the ice-rule cases compared with the random cases. For 
example, the van-Hove singularities and the semimctallic 
dip are retained up to fairly large U, while such structures 
are rapidly smeared out for the random case. 

It is noteworthy that these universal features of DOS 
are well captured by the exact solutions for the Husimi 
cactus models. The cactus models also show relatively 
small Uc for the transition to charge ice insulator; Uc = 2 
{Uc = 1) compared to the bare bandwidth W = 8 
{W = 6) for the tetrahedron (triangle) Husimi cactus. 
Moreover, at U oo, the cactus models show qualita- 
tively similar DOS as the original lattice models; a ID 
tight-binding form (plus a S function) for the tetrahe- 
dron (triangle) Husimi cactus. The agreement not only 
indicates that the cactus models give good references to 
the original lattice models but also confirms that these 
features of DOS are universal among the ice- rule systems. 

In contrast to DOS, transport properties are depen- 
dent on the lattice structures. For the pyrochlore lattice, 
the optical conductivity 17(0;), the low-energy weight w, 
and IPR consistently show that the system becomes insu- 
lating at U = Uc- This suggests that the metal-insulator 
transition occurs simultaneously with the transition to 
the charge ice insulator, and that it is not driven by 
the Anderson localization before the gap opening. In- 
terestingly, the low-energy weight indicates a nonmono- 
tonic change with increasing [/, presumably owing to the 
semimctallic nature at U = [Fig. HJb)] . 

The results for the checkerboard lattice look different. 
While the energy gap opens at Uc = 2.7(2) for this model, 
IPR indicates that the system is insulating already for 
U < 2, as shown in Fig. [TT] This appears to contra- 
dict with a finite w remaining up to U c:^ Uc', however, it 
is hard to estimate the true low-energy contribution be- 
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cause of the sharp dip structure of (t(cj) at cj ^ [Fig. [9]. 
Thus the results indicate that the Anderson insulator re- 
gion appears in a wide region of U before the energy gap 
opens at Uc ^ 2.7(2). The difference of transport prop- 
erties between the pyrochlore and checkerboard models 
may be attributed to the difference in dimensionality of 
the system. It is well known that the Anderson localiza- 
tion is relevant for lower-dimensional systems. 

From the viewpoint of dimensionality, the kagome 
model would share qualitative aspects of the transport 
properties with the checkerboard model. Our numerical 
results, however, do not clearly show the Anderson lo- 
cahzation (Fig. [Tn]). A possible reason for the apparent 
absence of the Anderson localization is that it occurs in 
a narrow range of U beyond our numerical resolution. 
The smaller Uc ^ 1-4(2) makes it harder to observe the 
Anderson transition if any. 

Recently, the Anderson localization in the systems 
with flat bands was analysed (ssj . It was shown that a 
random potential broadens the flat bands, and the states 
originating from the flat bands become critical, neither 
Anderson localized nor spatially extended. The analysis 
was done for, typically, a Gaussian distribution of the 
potential, which smears out the flat band divergence. In 
our case, the potential distribution is binary originat- 
ing in the configuration of localized particles, and hence, 
the divergences remain for finite amplitude of the po- 
tential. The problems are related with each other, but 
the situations are different. Our work focuses on the 
metal-insulator transition in the systems under the ice- 
rule correlated potentials; the results cast a new issue on 
the localization problems from the characteristic spatial 
correlation of potential configurations. 



V. SUMMARY 

In summary, we have studied the extended Falicov- 
Kimball model on the pyrochlore, checkerboard, and 
kagome lattices. With exact diagonalization of Hamil- 
tonian and approximation of the statistical average by 
the arithmetic mean, we obtained the density of states, 
optical conductivity, and IPR of the models. Through 
the analysis, we have clarified how the local ice-rule con- 
straint affects the global electronic structure, refiecting 
the peculiar nature of ice-rule manifold. 

The ice-rule local constraint gives rise to several uni- 
versal features in the density of states. The distinctive 
features are summarized as follows, (i) The energy gap 
opens at a much smaller U compared to the bandwidth, 
(ii) In the large U limit, the density of states becomes 
a similar form to that of one-dimensional tight-binding 
model (or a S functional peak), (iii) A cusp-like struc- 
ture appears as a precursor of the band-edge divergence 
in the large U limit. These features emerge from the in- 
teraction with localized particles under the ice rule, and 
are not observed in the coupling to random potentials. 
Furthermore, they are insensitive to the details of lattice 



structure, and can be captured by the exactly-solvable 
cactus models. The universal behaviors in the character- 
istic evolution of the density of states can be considered 
as a hallmark of transition to the charge-ice insulating 
state. 

In contrast, transport properties depend on the de- 
tailed lattice structure, such as dimensionality and lo- 
cal geometrical unit. For the pyrochlore lattice model, 
the optical conductivity and the inverse participation ra- 
tio consistently indicate that a metal-insulator transition 
takes place accompanied with the gap opening, in sharp 
contrast to the Anderson localization for the random 
case. The low-energy weight of the optical conductiv- 
ity shows a nonmonotonic behavior as a function of U 
in the metallic region. On the contrary, transport prop- 
erties for the checkerboard lattice model look different. 
Contrastive to the pyrochlore case, the low-energy weight 
decreases monotonically with respect to U. Moreover, 
the results of the inverse participation ratio and energy 
gap suggest Anderson localization in the region U < Uc- 
The contrasting behaviors are ascribed to the difference 
of dimensions of the systems. Yet different behaviors ap- 
pear for the kagome lattice models which have a different 
geometrical unit, a triangle: our results show no clear ev- 
idence of the Anderson localization for U < Uc- This is 
presumably due to the small Uc or the limited precision 
of our calculations. 

Our results clearly indicate the significance of local 
correlations on the global electronic structure and trans- 
port properties. Such local correlations should be re- 
sponsible for the peculiar pro perties observed in several 
pyrochlore compounds |13l4l5l |. where itinerant electrons 
interact with localized moments under the ice-rule type 
constraint. For understanding these interesting proper- 
ties, it will be crucial to take into account the charac- 
teristic spatial correlations emerging from the local con- 
straint. 
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Appendix A: Lifting of the ground state degeneracy 

In the main text, in replacing the statistical average 
by the arithmetic mean, we assumed that the Boltz- 
mann weight for all teh ice-rule configuration is equiv- 
alent. However, in the extended Falicov-Kimball model, 
Eq. IH) , the coupling to itinerant fermions may lift the de- 
generacy of ice-rule manifold and select a unique ground 
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state or a submanifold. Hence, to properly investigate 
the thermal properties of this model, we need to take 
account of the energy levels in the Boltzmann weight in 
Eq. ([S]). Nevertheless, we expect our results to be plau- 
sible at low temperatures as we will show below in this 
section. This is due to the fact that lifting of degeneracy 
is very small in the entire range of U, specifically, in the 
order of lO^'^t. 

To investigate the effect of coupling to itinerant 
fcrmions, we first consider the problem by the pertur- 
bation in t/U in the large U limit. For the pyrochlore 
lattice model, the lowest-order contribution to the en- 
ergy difference comes from the ring-exchange-type hop- 
ping process on a hexagon embedded in the pyrochlore 
lattice structure (see Fig. [T7|) . The contribution is in the 
sixth-order of t/U, and depends on the number of local- 
ized particles in the hexagon. The energy becomes the 
lowest for a hexagon with three localized particles. The 
relative energy to a fictitious reference system in which 
all hexagons have three localized particles is given by 

^E^JjE {12(^0 + ^e) + 10(ni + ng) + 20(^2 + rn)} , 

(Al) 

where nm (0 < m < 6) gives the number of hexagons 
in which m localized particles exist. Notes that the en- 
ergy depends only on the number of localized particules 
in each exagon, not on their configurations within each 
hexagon. Since the contribution starts from the order of 
t^/C/^, the energy difference among different ice- rule con- 
figurations remains very small in the large-?/ charge-ice 
insulating regime. 

To see how the energy difference develops as U de- 
creases beyond the perturbation regime, we numerically 
evaluate the energy difference in the entire range of U . 
Since the energy difference is very small and becomes 
comparable to the numerical resolution of our calcula- 
tions when we perform the arithmetic mean over the ran- 
domly generated samples, we try to estimate a typical 
difference by considering two specific periodic configura- 
tions with very different number of hexagons. One 
is the A-type stripe order shown in Fig. I17f a). in which 
localized particles exist only on [110] chains. The other 
is the hexagonal ordered state shown in Fig. [TTT b). in 
which one-third of hexagons in all the [111] kagome lay- 
ers are fully occupied by localized particles. These two 
states are largely different in the numbers of different- 
type hexagons as shown in TableU thus, the ground state 
energy is expected to be substantially different in the 
scheme of perturbation theory. We compare the ground- 
state energies of itinerant fermions for these two systems 
by numerical diagonalization method. The calculation 
was done on 16^ superlattices of 4 x 3'^ sites. 

The results are shown in Table [III From the numeri- 
cal results in a broader range of U including the metallic 
region, we find that the energy variance remains surpris- 



ingly small, typically in the order of 10 or less. 

The weak lifting of the ice-rule manifold implies that 




FIG. 17. (color online). Schematic pictures of (a) A-type 
stripe and (b) hexagonal ordered states. The sites with 
filled (open) circles represent —U/2 [-\-U/2) sites. The high- 
lighted bonds in (a) shows stripes of — [7/2 and (b) represents 
hexagons of —U/2 potential sites. See the text for details. 
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TABLE I. Distributions of hexagons with m localized par- 
ticles for A-type stripe and hexagonal ordered states. The 
values are normalized to the number of each hexagon in a 
cubic unit cell shown in Fig. [171 



the manifold is preserved down to very low T even in 
the ex-FK model: The Boltzmann weights for different 
ice-rule configurations are virtually the same for T /t ^ 
10"'^. On the other hand, the ice-rule constraint of the 
localized particles is expected to be well retained for T < 
U . In such wide T region, thermal average might be well 
approximated by the arithmetic mean as performed in 
the present study. Such situation was considered also in 
our previous study [23 |. 

The situation is similar for the kagome lattice models. 
In this case also, the smallest loop is the six-site hexagon, 
and hence, the discussions above can be applied straight- 
forwardly. 

For the checkerboard lattice, the same approach gives 
fourth-order perturbation to be the lowest order 

A£;= — {2(no + ?i4) + 6n2}, (A2) 

since the smallest loop is on the four-site plaquette [see 
Fig. lUb)]. This gives larger energy shift compared to 
the pyrochlore and kagome lattice cases, and hence, is 
expected to be more relevant on lifting the ground state 
degeneracy. For this reason, some ordcrings may be ex- 
pected for the checkerboard lattice models at low T. Nev- 
ertheless, in this paper, we limit ourselves to treating this 
model with the arithmetic mean and leave the possibility 
of ordering for future problem. 
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-3.20546 


-3.20474 


-7.2 X 10"" 
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-2.67592 


-2.67658 


6.7 X 10"'' 


10 


-2.37914 


-2.37924 


0.9 X IQ-" 



TABLE II. Ground state energy for the A-type stripe and 
hexagonal ordered states per site and their difference. The 
data are for 16'' superlattices of 4 x 3^ sites. 
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